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A method to produce a spin-dependent rotating potential using near-resonant circularly polarized 
laser beams is proposed. It is shown that half-quantum vortices are nucleated in a spinor Bose- 
Einstein condensate with an antiferromagnetic interaction. In contrast to the vortex nucleation in a 
scalar BEC, the spin- vortex nucleation occurs at a low rotation frequency (— O.Icj^ with oj± being 
the radial trap frequency) in a short nucleation time ~ 50 ms without dissipation. A method for 
nondestructive measurement of half-quantum vortices is discussed. 
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I. INTRODUCTION 

Quantized vortices are hallmarks of superfluidity and 
have widely been studied in Bose-Einstein condensates 
(BECs) of atomic gas. In scalar BECs, vortex states 
have been generated using a variety of methods, e.g., 
phase imprinting p], 0, 0] and potential rotation 0, Q . 
Vortices in BECs of atoms with spin degrees of freedom 
(spinor BECs) have been attracting increasing interest 
since the Berkeley group [6] recently observed sponta- 
neous creation of polar-core vortices in the magnetiza- 
tion of a spin-1 87 Rb BEC using a nondestructive spin- 
sensitive imaging technique [7]. In this experiment, spin 
vortices are created through the dynamical instability 
arising from ferromagnetic interaction The MIT 

group has generated Mermin-Ho and Anderson- Toulouse 
vortices in a spin-1 23 Na BEC using the Berry phase- 
imprinting method [9] . These spin- vortex states are pre- 
dicted to be stable in a rotating potential [l0[. However, 
spin- vortex nucleation by a rotating potential has not yet 
been realized experimentally. 

For a scalar BEC in a rotating potential, there is a crit- 
ical rotation frequency for the vortex nucleation [4|, , 
above which the surface mode becomes dynamically un- 
stable [l2|, Qjjj], allowing vortices to enter the condensate. 
According to the numerical simulations in Ref. P3], en- 
ergy dissipation is needed to reproduce the dynamics of 
vortex nucleation observed in the experiments. 

In the present paper, we propose a method to create 
spin vortices using a spin- dependent rotating potential, 
which is produced by near-resonant circularly polarized 
laser beams. Using this method, we can rotate spin sub- 
levels selectively. We numerically demonstrate that half- 
quantum spin vortices [HI, EH Q3 are nucleated in the 
antiferromagnetic ground state of the spin-1 23 Na BEC. 
The critical rotation frequency for this spin- vortex nucle- 
ation is O.la;^ with co± being the radial trap frequency, 
and the nucleation occurs at t ~ 50 ms without dissipa- 
tion. This is in marked contrast to the vortex nucleation 
in a scalar BEC, in which the critical rotation frequency 
is ~ 0.7cj [IH, Qjjj] and vortices are never nucleated for 
t < 100 ms unless dissipation is taken into account [lij . 



A method to observe half-quantum vortices in a nonde- 
structive manner is also proposed. 

This paper is organized as follows. Section HJ provides 
a method to create a spin-dependent optical potential. 
Section IIIII numerically demonstrates the nucleation of 
half-quantum vortices in a spin-1 23 Na BEC. Section HVl 
discusses a method for nondestructive measurement of 
half-quantum vortices. Section [V] presents our conclu- 
sions. 



II. SPIN-DEPENDENT OPTICAL POTENTIAL 

It is known that a far-off-resonant laser beam produces 
a potential that is independent of for the hyperfine 
state \F : tuf) [18]- We show here that an appropriately 
tuned laser beam can produce a potential that strongly 
depends on m^. For simplicity, we restrict ourselves to 
alkali atoms with electron spin 5=1/2 and nuclear spin 
1 = 3/2 (e.g., 23 Naand 87 Rb). 

We consider an atom in the electronic ground state 
|g, F = l,ra F ) located in circularly polarized laser field 
with frequency c^o- The energy shift due to the ac Stark 
effect for the ground state is proportional to 



|(n|(7±|g,F = l,m F ) 
Hluq — E n 



(i) 



where the summation is taken for the relevant states \n) 
with energy E n . The dipole operators a± in Eq. (pQ) 
change the orbital angular momentum of the outermost 
electron from \L = 0, tul = 0) to \L = 1, tul = ±1). 

When the laser frequency is close to the Di and D 2 
lines, Eq. (pQ) is approximated to be 
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D2 transitions, the lifetime of the BEC is shortened by 
the spontaneous emission. We assume that the domi- 
nant contribution to the trapping potential is made by 
the far-off-resonant beam (with intensity ex Af ar and de- 
tuning Af ar ), and the a± polarized beams (with intensity 
oc A near and detuning A near ) are only small perturbations 
[say, (A near /A n ear )/(-4far /Af ar ) = 0.05]. The ratio of the 
loss rates then becomes 



r /A 2 



^far/Af ar 



0.05- 



ifar 



(4) 



In Ref. [IK, Afar ~2x 10 14 Hz, and the loss rate for 
23 Na atoms is c± 0.03 Hz. Using these parameters, the 
right-hand side of Eq. (|4|) is ~ 40, and the lifetime for 
the present system is estimated to be ~ 0.8 s. 



III. DYNAMICS OF SPIN VORTEX 
NUCLEATION 

A. Formulation of the problem 



FIG. 1: Transition strengths with circularly polarized laser 
beams for the (a) Di and (b) D2 lines, (c) Subtraction of the 
transition strengths in (b) from those in (a) for the case of 
hevo = (Ed, +£d 2 )/2. 



where E^ ± and E^ 2 are the energies and |Di, F f , mp>) 
and |D2, F\ mp/) are the excited states for the Di and 
D2 lines, respectively. Since J = 1/2 (J = 3/2) and 
/ = 3/2 for the Di (D2) state, the possible hyperfine 
spins are F f = 1,2 (F' = 0,1,2,3), where J = L + S. 
Calculating each term in Eq. ([2]), we obtain the transition 
strengths given in Figs. [1] (a) and[T] (b). 

If, in particular, the laser frequency is tuned to the 
center of the Di and D2 lines, i.e., Hluq = (Ed 1 +£?d 2 )/2, 
the factors before the summations in Eq. ([2]) have the 
same magnitude with opposite signs. The transition coef- 
ficients in this case are shown in Fig.[T](c). We note that 
the transition coefficients for mp = ±1 — > mp> = ±2 
vanish. This indicates that a cr + (cr_) laser does not af- 
fect the mp = 1 (—1) state. The produced potential is 
an attractive potential, since the coefficients are negative. 
Thus, the <j± polarized beams with hcoo = (Ed 1 +^d 2 )/2 
combined with a far-off-resonant beam produce an opti- 
cal potential V mF (r) for each mp state as 



V(r) + -[V+(r)+V-(r)]. 



V (r) 

V ±1 (r) = V(r) + V T (r), 



(3a) 
(3b) 



where V(r) and V±(r) are proportional to the strengths 
of the far-off-resonant beam and the a± polarized beams, 
respectively. 

We estimate the lifetime of the BEC. Since the laser 
frequency ujq = (Er> ± + E^ 2 )/(2h) is close to the Di and 



We employ the zero-temperature mean-field approxi- 
mation. The dynamics of a BEC for spin-1 atoms are 
described by the three-component Gross-Pitaevskii (GP) 
equations given by 



ih 



dt 



h 

— V 2 + Fo(r,t)+c n 



^0 



-^(F + ^i+F-^-i), 



(5a) 



ih, 



dt 



2M 



V 2 + V ±1 (r,t)+c n 



+ci I ~^ F T^o ± F x il>±i 



±1 



(5b) 



where M is the mass of an atom, and ip m describes the 
mean-field wave functions satisfying 



]T j \<p m \ 2 dr = N 



(6) 



with N being the number of atoms. The interaction co- 
efficients cq and c\ are given by 



Airh do + 2a2 
~M 3 



Cl = ~M — ' (7) 



where ao and a 2 are the s-wave scattering lengths for 
colliding channels with total spins and 2, respectively. 
The system is ferromagnetic for c\ < 0, and antiferro- 
magnetic or polar for c\ > 0. In Eq. ([5]), the atomic 
density n is defined as 



(8) 



m=— 1 
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and the magnetization has the forms 

F z = |Vi| 2 -|V-i| 2 , (9) 
F + = F*_ = >/2 (ViVto + V'oV'-i) • (10) 

In the present analysis, we consider 23 Na atoms, which 
have an antiferromagnetic interaction 20]. We use the 
scattering lengths ao + 2<22 = 53.4<2b [2l| and a2 — «o = 
2.47aB [22|, where an is the Bohr radius. For the ini- 
tial state, we first prepare the antiferromagnetic ground 
state ip±i = V'ini and = using the imaginary-time 
propagation method. If = in the initial state, the 
right-hand side of Eq. ([5a|) vanishes, and always re- 
mains 0. In realistic situations, however, quantum and 
thermal fluctuations and residual atoms may trigger the 
growth of the m = component. We therefore simulate 
this possibility by giving the small white noise to the 
initial state of ipo as 

^oW=ATe(r), fe(r)=^ini(r), (11) 

where AT is a normalization constant and e(r) includes 
the complex random numbers obeying the normal distri- 
bution e~' e ' /( 2cr ^ /(2ira 2 ). The random numbers are set 
to each point of the numerical mesh. The value a is cho- 
sen to be (7 — 3.5 x 10 -3 , for which the initial population 
of the m = component is about 1 %. 

In the initial-state preparation, only the far-off- 
resonant laser beam is applied, which produces the spin- 
independent axisymmetric trapping potential, 

V(r) = ^M[co 2 ± (x 2 +y 2 )+co 2 z 2 }. (12) 

For t > 0, we additionally apply a time-dependent cr_ 
polarized laser beam, producing a potential rotating at a 
frequency ^, 

V.(r,t) = ~Mw' 2 X(t) 2 , (13) 

where 

X(t) = xcosQt + y sin Qt. (14) 

This potential can be generated using two beams rotating 
around the center of the trap. In the following calcula- 
tion, we take uo' 2 = 0.050^. We do not apply the <r+ 
polarized laser beam, i.e., 

V+(r)=0. (15) 

The m = — 1 component therefore does not undergo a 
rotating potential [see Eq. (|3J)]. 

We assume that the trapping potential V(r) is tight in 
the z direction and the system is effectively two dimen- 
sional (2D). When %uj z is much larger than the charac- 
teristic energy of the system, the wave function in the z 
direction is frozen in the ground state of the harmonic po- 
tential. The effective 2D interaction strength is obtained 
by integrating the GP energy functional with respect to z. 



We use the trap frequencies (u;±,(j z ) = 2tt x (120,5000) 
Hz. 

The time evolution of the system is obtained by nu- 
merically solving the 2D GP equation using the finite 
difference method with the Crank-Nicolson scheme. We 
divide 38.2 fim x 38.2 /mi space into a 200 x 200 mesh. 
We have verified that the results do not depend on the 
size of the mesh. 



B. Nucleation of half-quantum vortices 

Figure [2] shows the time evolution of the density and 
phase profiles with TV = 2 x 10 4 . The rotation fre- 
quency of the potential V- in Eq. ([T3j) is chosen to be 
Q = 0.13c<;^. At t = 20 ms, both rn = ±1 components 
start to deform, and the two topological defects approach 
the m = 1 component. At t = 50 ms, the two topological 
defects enter the m = 1 cloud, where the density holes 
in the rn = 1 component are filled with the rn = — 1 
component. After that, the topological defects leave the 
condensate (t = 150ms). Interestingly, the topological 
defects enter the condensate again (t = 250 ms), and 
the entry-exit cycles are repeated. The total density is 
almost unchanged throughout the dynamics because of 
Co ^> c\. We find that no appreciable spin-exchange dy- 
namics occurs, and the m = component remains small 
(< 3 %) for t < 300 ms. 

The topological spin structures in Fig. [2] (at t = 50 ms 
and 250 ms) are the half-quantum vortices [HI, [TH, 03 • 
In a spin-1 system, the general form of the half-quantum 
vortex located at x = y = is given by [l6[ 

/ V'l \ //i(r±)e ± ^\ 
*hqv = V>o = , (16) 

VV-i/ V /-i(rx) J 

where r± = (x 2 + y 2 ) 1 ^ 2 and (j) = arg(x + iy). The func- 
tion /i(rj_) vanishes at r± = and in a infinite system 
/i(oo) = /_i(oo) should be satisfied. The spatial ro- 
tation of ^hqv around the z axis is related to the spin 
rotation as 

e-^* hqv = e^/ 2 e^*/ 2 * hqV) (17) 

where L z = —id^, F z = m, and x is an arbitrary angle. 
Equation ([TT]) indicates an interesting fact: when we go 
around the half-quantum vortex core (x = 27t), the spin 
state rotates only by ±7r. Thus, the half-quantum vortex 
has a structure similar to a Mobius strip. The spin vor- 
tices shown in Fig. [2] have the same topological structure 
as Eq. ([T6]) in the vicinity of the vortex cores. We can 
therefore regard the dynamics in Fig. [2] as half-quantum 
vortex nucleation. 

For a scalar BEC, the critical rotation frequency above 
which vortex nucleation occurs is ft ~ 0.7uu_\_, and the 
typical nucleation time is ~ 100 ms [11]. In order to 
reproduce this nucleation time by the GP equation, we 
must take into account the effect of dissipation [l4|. In 
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FIG. 2: (Color) Time evolution of the density and phase pro- 
files of the m = ±1 components and the total density profile 
for the potential given in Eqs. ([T2]) - ([T5]) with f2 = 0.13u±. 
The m = 1 component is affected by the rotating potential. 
The m — component is negligibly small. The field of view 
of each panel is 38.2 /mi x 38.2 /mi. The unit of the density 
is N/al Q with a ho = [h/(Mcu ± )} 1/2 . 

contrast, the spin- vortex nucleation in the present system 
occurs at much lower rotation frequency Q = 0.13uj_\_, and 
the nucleation time is t ~ 50 ms even without dissipa- 
tion. These significant differences between the scalar and 
spinor systems originate from the energy cost for vortex 
nucleation. For a scalar BEC, the energy cost is deter- 
mined by Co because of the density hole at the vortex 
core, whereas the energy cost by the core of the spin 
vortex is determined by c\ (< Co). 

Figure [3] shows time evolution of the orbital angular 
momentum in the m = 1 component, 




We note that L\ remains small for Q = 0.05co± and 
0.5o;^, whereas L\ becomes large at Vt = 0.13o;^. This 
implies that there is a region in which the dynamical 
instability against spin-vortex nucleation sets in. The 
oscillation of the green curve in Fig. [3] corresponds to 
the cycles of entry and exit of the half-quantum vortices 
shown in Fig. [2j 

For ft > 0.2cjx, we found that the vortices nucleate not 
only in the m = 1 component but also in the m = — 1 
component even though V-i is not a rotating potential. 
This is because the m = — 1 component effectively feels 
a rotating potential through the interaction with the ro- 
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FIG. 3: (Color) Time evolution of the orbital angular momen- 
tum Li in the m — 1 component for Q = 0.05o;^, 0.13o;^, and 
0.5o;^. The insets show the density and phase profiles at 50 
ms for Q = 0.13o;± and at 20 ms for Q = 0.05uu±. 

tating m = 1 component. By the same mechanism, we 
can create the vortices only in the m = — 1 component 
when V- in Eq. ([T3]) is positive. The dependence of the 
dynamics of the spin-vortex nucleation on the external 
potential merits further study. 

We have considered so far the case of zero magnetic 
field. When the magnetic field is applied in the z direc- 
tion, the m = component is energetically favored and 
its growth is enhanced. In order to suppress the growth 
of the m = component, the magnetic field must be 
B <J 10 mG. When the magnetic field is B = 10 mG and 
the initial population of the m = component is ^ 1 
%, the m = population is suppressed below 3 % for 
t < 300 ms. 



IV. NONDESTRUCTIVE MEASUREMENT OF 
HALF-QUANTUM VORTICES 

Using the nondestructive spin-sensitive imaging tech- 
nique developed by the Berkeley group [3] , we can observe 
the half-quantum vortices. In the Berkeley method, a <r + 
circularly polarized probe light is shone in the y direc- 
tion and phase-contrast images are obtained. The phase- 

(2) 

contrast signal is proportional to Aon + A\F y + A^Fy , 
where F y = (F+ + F-)/(2i) [see Eq. CEOj)] and 

F v 2) = E C(/^)wVw = |^i-^-i| 2 /2+|^o| 2 (19) 

771,777/ 

with f y being the y component of the spin-1 matrix, and 
the coefficients Aq, A\, and A^ are given in Tab. H In 
the Berkeley experiments @, 0, the F = 1 -> F' = 2 Di 
transition was used. For a ferromagnetic BEC, the mag- 
nitude and phase arg(F + ) of the transverse magne- 
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TABLE I: Coefficients of the phase- contrast signal oc Aon + 

( 2) 

Al-Fj/ + ^Fy for each transition. Signal-to-bias ratio is pro- 



portional to R = A 2 /(A + A 2 /2). 



tization are obtained from the oscillation of F y due to 
the Larmor precession. 

In the present case, the spin state is written as 



(Ci(r)e- 


VC-i( r )' 



(20) 



where cjl is the Larmor frequency and £±i depends only 
on r in the time scale of ou^ . For this state, F y ~ 0, and 
the phase-contrast signal is proportional to 

A n + \A 2 [|Ci | 2 + |C-i I 2 + 2Re (CiC-ie 2 ^)] ■ (21) 

From this signal oscillating at the frequency 2cjl, we 
can determine the spatial distribution of the relative 
phase between the m = ±1 components. Around a 
half-quantum vortex, the phase of the oscillating signal 
changes by 2tt. The ratio of the oscillating signal to the 
bias is proportional to 



R 



A,, 



(22) 



A +A 2 /2 

For the present purpose, therefore, the F = 1 — > F' = 



D2 transition may be most suitable, since R is largest 
(Tab.[I|). 



V. CONCLUSIONS 

We have proposed a method to create a spin-dependent 
optical potential using near-resonant circularly polarized 
laser beams. We have shown that spin vortices can be 
nucleated in a spinor BEC using the spin-dependent ro- 
tating potential. 

We considered a situation in which only the m = 1 
component of the antiferromagnetic ground state of a 
spin-1 23 Na BEC is stirred by the spin-dependent rotat- 
ing potential, and found that half-quantum vortices enter 
the condensate. To our knowledge this is the first pro- 
posal for nucleating fractional vortices by a rotating stir- 
rer. The spin- vortex nucleation occurs at a low rotation 
frequency (~ 0.1o;_l), compared with vortex nucleation 
in a scalar BEC (~ 0.7lj±). Moreover, the spin vortices 
easily enter the condensate: the nucleation time is t c± 50 
ms without dissipation. The spin vortices exit from the 
condensate, and the entry-exit cycles are repeated. We 
have also shown that the half-quantum vortices can be 
observed in a nondestructive manner using the method 
of the Berkeley group. 

The spin-dependent optical potential is a powerful tool 
for manipulating a spinor BEC, and may be applied to 
the generation of various spin textures. 
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